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(57) Abstract: The invention relates to a device and a method for analysing a sample plate on which at least two material samples 
are arranged. According to said method, an impedance spectrum is measured for each material sample, and a structure of a switching 
circuit equivalent comprising at least one electronic component is then determined according to the respectively measured impedance 
spectrum. Starting values for the components of the respective switching circuit equivalent are determined for an error minimising 
calculation. During the error minimising calculation, a theoretical impedance spectrum for at least one of the material samples is 
calculated on the basis of the impedance spectrum measured for the material sample, the starting values for the components of the 
^^^1 corresponding switching circuit equivalent, and fit values for the components of the corresponding switching circuit equivalent. A 
^ validation variable is then determined for the calculated theoretical impedance spectrum, and an evaluation variable is determined 
by comparing at least one of the fit values for the components with a reference value. 

[Fortsetzung aufder ndchsten Seite] 
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(57) Zusammenfassung: Es werden eine Vorrichtung und Verfahren zur Analyse einer Probenplatte vorgeschlagen, auf der mindes- 
tens zwei Materialproben angeordnet sind. Bei dem Verfahren wird fiir jede der Materialproben ein Impedanzspektrum gemessen. 
In Abhangigkeit von dem jeweils gemessenen Impedanzspektrum wird ein Aufbau eines Schaltkreisaquivalents bestimmt, der min- 
destens ein elektronisches Bauelement umfasst. Fiir eine Fehlerminimierungsrechnung werden dann Startwerte fiir die Bauelemente 
des jeweiligen Schaltkreisaquivalents bestimmt- Bei der Fehlerminimierungsrechnung wird ein theoretisches Impedanzspektrum 
fiir mindestens eine der Materialproben unter Zugrundelegung des fiir die Materialprobe gemessenen Impedanzspektrums sowie der 
Startwerte fiir die Bauelemente des betreffenden Schaltkreisaquivalents berechnet und Fitwerte fiir die Bauelemente des betreffenden 
Schaltkreisaquivalents bestimmt. Anschliessend wird eine Validierungsgrosse fiir das errechnete, theoretische Impedanzspektrum 
bestimmt und eine Auswertegrosse durch Vergleich mindestens eines der Fitwerte fiir die Bauelemente mit einem Referenzwert er- 
mittelt. 
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Vorrichtung und Verfahren zur Analyse einer Probenplatte 
Stand der Technik 

Die Erfindung betrifft eine Vorrichtung zur Analyse einer 
Probenplatte gemali der im Oberbegriff des Patentanspruches 
1 naher definierten Art sowie ein Verfahren zur Analyse ei- 
ner Probenplatte. 

Insbesondere im Bereich der Materialwissenschaf ten, der 
Chemie und der Pharmazie ist es von erheblicher Bedeutung, 
hinsichtlich des jeweiligen Anwendungsf alls optimierte 
Stoffe und Materialien aufzufinden bzw. zu entwickeln, Ein 
spezielles Anwendungsgebiet stellt hierbei die Sensorik 
dar, die eine wichtige Schlusseltechnologie mit einer ste- 
tig wachsenden Zahl von Anwendungen sowohl im indust riellen 
als auch im privaten Bereich ist. So werden Sensoren bei- 
spielsweise bei technischen Prozessuberwachungen, im Be- 
reich des Umweltschutzes, im Bereich der Medizin oder auch 
im Kraf tf ahrzeugbereich eingesetzt. Ein betracht liches Maft 
an Entwicklungsarbeit wird derzeit insbesondere in die Ent- 
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wicklung schneller und hochsensitiver Sensoren mit einer 
geringen Querempf indlichkeit gesteckt . 

Die Entwicklungen beschrankten sich bisher in der Regel auf 
eine Optimierung bzw. Modif izierung bekannter Materialien. 
Es besteht jedoch das Problem, dass fur bestimmte Anwen- 
dungsbereiche der Sensorik ein hoher Bedarf an neuen Mate- 
rialien besteht, der mit konventionellen Methoden, die sich 
durch eine Herstellung von Einzelsensoren und eine an- 
schlieBende sequent ielle Charakterisierung auszeichnen, 
nicht hinreichend gedeckt werden kann. 

Insbesondere bei der Entwicklung neuartiger, sensitiver Ma- 
terialien bzw. Materialkombinat ionen kann es zweckmaBig 
sein, Verfahren aus dem Bereich der kombinatorischen Chemie 
bzw. sogenannte High-Throughput-Methoden einzusetzen, Bei 
diesen Verfahren handelt es sich um parallelisierte Synthe- 
se- und Screeningverf ahren, durch welche sich neue Materia- 
lien bzw. Materialkombinationen erschliel3en lassen konnen 
bzw, bereits bekannte Syntheseverf ahren fiir existierende 
Materialien in einem breiten Parameterf eld optimiert werden 
konnen . 

Aus der US 5,985,356 ist eine allgemeine Darstellung von 
High-Throughput-Verf ahren bekannt, wobei in dieser Schrift 
insbesondere die Anwendung der im Wesentlichen aus dem Be- 
reich der Pharmazie bekannten kombinatorischen Chemie auf 
chemische und ma terialwissenschaf tliche Anwendungsgebiete 
vorgeschlagen wird . 

Eine Vorrichtung der einleitend genannten Art zur Analyse 
einer Probenplatte ist beispielsweise aus der 
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DE 101 31 581 Al bekannt. Diese Vorrichtung umfasst eine 
Probenplatte, auf der 64 Materialproben matrixartig aufge- 
bracht sind, welche jeweils mit zwei Elektroden verbunden 
sind^ die wiederum mit Kontaktstellen versehen sind, an 
welchen ein Mittel zur reversiblen und adressierbaren Kon- 
taktierung zur Anlage gebracht werden kann. 

Vorteile der Erfindung 

Die erf indungsgemaiie Vorrichtung zur Analyse einer Proben- 
platte mit den Merkmalen nach dem Oberbegriff des Patentan- 
spruches 1 und mit einem in einen Gehausetrager einsetzba- 
ren Messkopf, der zur elektrischen Verbindung mit den Kon- 
taktiermitteln je Materialprobe zwei Messdrahte umfasst, 
die mit Vorspannung an Kontaktf lachen der Probenplatte ' an- 
liegen und mit einer Mess- und Auswerteeinheit verbunden 
sind, hat den Vorteil einer einfachen Handhabung, da bei 
Verwendung einer genormten Probenplatte der Kontakt zwi- 
schen den Materialproben und der Mess- und Auswerteeinheit 
durch einfaches Einsetzen des Messkopfes in den Gehausetra- 
ger herstellbar ist. 

Die Vorrichtung nach der Erfindung eignet sich insbesondere 
zur Entwicklung und Auffindung von Materialien bzw. Materi- 
alkombinationen, die als Sensormaterialien einsetzbar sind 
und uber ihre elektrischen Eigenschaf ten charakterisierbar 
sind. Beispielsweise kann die Vorrichtung nach der Erfin- 
dung zur Entwicklung eines optimierten Sensormaterials ei- 
nes Gassensors eingesetzt werden. 
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Mit der Vorrichtung nach der Erfindung ist es moglich, eine 
hohe Anzahl an potentiellen Sensormaterialien^ die auf der 
Probenplatte angeordnet sind, unter verschiedenen Prufgasen 
bei unterschiedlichen Temperaturen, die beispielsweise bis 
zu 800 °C betragen konnen, nahezu simultan zu untersuchen. 
Die Untersuchung kann potentiometrisch, resistiv^ kapazitiv 
Oder auch uber die komplexe Impedanzspektroskopie erfolgen. 

Der Begriff Spektroskopie ist im vorliegenden Fall so zu 
verstehen, dass es sich urn f requenzabhangige Messungen han- 
delt, d- h. dass die Impedanz einer Probe bei unterschied- 
lichen Messf requenzen untersucht wird. Beispielsweise wer- 
den die einzelnen Materialproben jeweils in einem Frequenz- 
bereich zwischen 10 Hz und lO"^ Hz mit einer Messdatendichte 
von 15 Messpunkte pro Dekade untersucht. Dies bedeutet, 
dass pro Mater ialprobe bei ' einer derartigen Messdatendichte 
180 Messdaten ermittelt werden. Eine weit reichende Datenre- 
duktion kann beispielsweise durch eine Anpassung eines ge- 
eigneten Schalt kreisaquivalents an die Messdaten erreicht 
werden . 

Urn die Messdrahte, die insbesondere uber Schmelzkugeln an 
den Kontaktf lachen der Probenplatte anliegen, unter der 
Vorspannung zu halten, konnen die Messdrahte jeweils mit 
einem insbesondere vergoldeten Federkontakt verbunden sein, 
der einen konstanten Anpressdruck des jeweiligen Messdrahts 
auf der jeweiligen Kontaktf lache gewahrleistet . 

Urn zur Entwicklung eines Gassensors die Materialproben mit 
einem bestimmten Prufgas beauf schlagen zu konnen, kann der 
Messkopf mit einer Gasversorgungseinheit verbunden sein. 
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Urn die Materialproben unterschiedlichen Prufgas- bzw. Refe- 
renzgasatmospharen aussetzen zu konnen, umfasst die Gasver- 
sorgungseinheit , die zweckma/iigerweise mit einer Datenver- 
arbeitungseinheit der Mess- und Auswerteeinheit verbunden 
ist, eine Gasmischeinrichtung . Ferner kann die Gasversor- 
gungseinheit zur Anfeuchtung der Prufgas- bzw. Referenzgas- 
atmosphare ein Wasserreservoir umfassen. 

Der Messkopf kann des Weiteren so ausgelegt sein, dass er 
als integrierten Bestandteil einen Gasraum umfasst, der 
uber den Materialproben der Probenplatte angeordnet ist und 
vorzugsweise von einer im Wesentlichen glockenf ormigen Ver- 
teilereinrichtung gebildet ist. Der Gasraum ist mit der 
Gasversorgungseinheit verbunden . 

Um eine homogene Verteilung des Pruf- bzw. Ref erenzgases in 
dem Gasraum zu erreichen, ist bei einer vorteilhaf ten Aus- 
fuhrungsform der Vorrichtung nach der Erfindung in dem Gas- 
raum ein Diffusor angeordnet. 

Um bei einer Vielzahl von Materialproben auf der Proben- 
platte die einzelnen Materialproben auf einfache Art und 
Weise vermessen zu konnen, umfasst die Mess- und Auswerte- 
einheit vorteilhafterweise zwei Relaisschaltf elder , die mit 
den Messdrahten verbunden sind und beispielsweise bei 64 
Materialproben auf der Probenplatte vorzugsweise jeweils 
eine 3x64-Matrix aus hochf requenztauglichen Relais aufwei- 
sen. In diesem Fall sind die 64 Materialproben auf der Pro- 
benplatte in einem Messzyklus vermessbar, wobei mindestens 
drei Messgroften zuganglich sind, und zwar beispielsweise 
uber einen Impedanz-Analysator die Impedanz und uber ent- 
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sprechende weitere Messgerate die Gleichstromwiderstande 
der Materialproben und deren Strom/Spannungs-Kennlinien . 

Die Mess- und Auswerteeinheit ist zweckmaliigerweise mit ei- 
ner Mess- und Steuersof tware ausgestattet , die einerseits 
den Messablauf steuert und andererseits gewonnene Messdaten 
an eine entsprechende Datei oder auch an eine relationale 
Datenbank iibergibt, die von einer Auswertesof tware eingele- 
sen werden kann. 

Die Auswertesof tware arbeitet vorzugsweise derart, dass sie 
eine Fitf unktionalitat zur Berechnung theoret ischer Impe- 
danzspektren fiir die einzelnen Proben umfasst^ wobei die 
Berechnung vorzugsweise unter Zugrundelegung eines Schalt- 
kreisaquivalents erfolgt, das mindestens ein virtuelles o- 
der reales elektronisches Bauelement umfasst. Ein virtuel- 
les Bauelement ist beispielsweise ein Constant Phase Ele- 
ment (CPE) . Die Fitf unktionalitat berechnet also auf Grund- 
lage der Messdaten fur ein beispielsweise aus einem seriel- 
len RC-Glied bestehendes Schalt kreisaquivalent ein an ein 
gemessenes Impedanzspektrum bestmoglich angenahertes , theo- 
retisches Impedanzspektrum, wobei zur Anpassung des theore- 
tischen an das gemessene Spektrum eine Variation der Kapa- 
zitat und/oder des Widerstands der Bauelemente des RC- 
Glieds durchgefuhrt wird. Dem bestmoglich angepassten, the- 
oretisch Impedanzspektrum sind also in diesem Fall ein Wi- 
derstandswert und ein Kapazitat swert fur das RC-Glied zuge- 
ordnet. Aus dem Wider standswert kann beispielsweise durch 
Vergleich mit einem Referenzwert auf die Sensitivitat der 
jeweiligen Mater ialprobe geschlossen werden. 
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Um bei einer Vielzahl von Materialproben^ die unter ver- 
schiedenen Messbedingungen vermessen warden, die fur die 
einzelnen Materialproben gewonnenen, beispielsweise die 
Sensitivitaten darstellenden Ausgabewerte einer einfachen 
Auswertung zuganglich zu machen, umfasst die Auswertesof t- 
ware vorteilhaf terweise eine Datamining-Funktionalitat . Die 
Datamining-Funktionalitat ermittelt beispielsweise unter 
Anwendung vorzugsweise mehrdimensionaler Zielf unktionen die 
fur den jeweiligen Anwendungsf all optimale Materialprobe 
auf numerischem Wege. 

Alternativ oder zusatzlich kann die Datamining-Funktionali- 
tat auch eine Visualisierungsf unkt ionalitat umfassen. In 
diesem Fall kann ein Benutzer bildschirmgestut zt die fur 
den jeweiligen Anwendungs f all optimale Materialprobe ermit- - 
teln. Beispielsweise arbeitet die Visualisierungsf unkt iona- 
litat mit einem Farbspektrum, wobei einer bestimmten Farbe 
eine hohe Sensitivitat der Materialprobe und einer anderen 
Farbe eine geringe Sensitivitat zugeordnet ist. 

Um auch Sensormaterialien auffinden zu konnen, die fur An- 
wendungsf alle geeignet sind, bei denen - wie beispielsweise 
bei einem Abgas - hohe Temperaturen herrschen, weist die 
Vorrichtung nach der Erfindung vorzugsweise eine Heizein- 
richtung auf, in welche die Probenplatte vorzugsweise ein- 
tauchbar ist. In diesem Fall weist die Vorrichtung insbe- 
sondere einen Hochtemperaturreaktor auf, der von der Heiz- 
einrichtung begrenzt ist und in dem die die Materialprobe 
umfassende Probenplatte mit unterschiedlichen Pruf- bzw. 
Ref erenzgasen beaufschlagt werden kann. 
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Die Erfindung hat auch ein Verfahren zur Analyse einer Pro- 
benplatte mit den Merkmalen gemali Anspruch 19 zum Gegen- 
stand. Durch Anwendung dieses Verfahrens ist as moglich^ 
eine Probenplatte, auf der eine grolie Anzahl an Material- 
proben, beispielsweise 64 Materialproben angeordnet sind, 
unter unterschiedlichen Bedingungen elektrisch zu vermessen 
und in vollautomatischer Weise eine fur den betreffenden 
Anwendungsf all am besten geeignete Materialprobe auszuwah- 
len. Dies erfolgt durch automatische Auswahl von Startwer- 
ten fur die Bauelemente des jeweiligen Schalt kreisaquiva- 
lents und die anschlieftende Fehlerminimierungsrechnung . Die 
Startwerte werden bei der Fehlerminimierungsrechnung einge- 
setzt, wobei ausgehend von den Startwerten ein an das je- 
weils gemessene Impedanzspektrum angepasstes, theoret isches 
Impedanzspekt rum fur die betreffende Materialprobe unter 
den jeweiligen Messbedingungen berechnet wird. Durch dieses 
Vorgehen ist eine weitgehende Datenredukt ion moglich, da 
ausgehend von den Impedanzspektren, die eine Vielzahl an 
Messpunkten haben, wenige abgeleitete Grofien darstellende 
Fitwerte ermittelt werden, die die einzelnen Bauelemente 
des Schaltkreisaquivalents beschreiben. Die Fitwerte stel- 
len also Dimensionierungen der Bauelemente des Schaltkrei- 
sesaquivalents dar, mit denen das gemessene Impedanzspekt- 
rum bestmoglich simuliert werden kann. 



Bei der Bestimmung des Schaltkreisaquivalents kann bei- 
spielsweise eine serielle Schaltung von vier RC-Gliedern 
ausgewahlt werden, wobei die hoheren RC-Glieder gegebenen- 
falls auf einen Startwert gesetzt werden, der keinen Ein- 
fluss auf die Berechnung des theoretischen Impedanzspekt- 
rums hat . 
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Die fur die Fehlerminimierungsrechnung er f orderlichen 
Startwerte eines RC-Glieds werden vorzugsweise aus der ma- 
ximal gemessenen, imaginaren Impedanz Z"_MAX und der kor- 
respondieren Messfrequenz f_Z"_MAX gema/i folgenden Formeln 
berechnet : 

R\_START= -2.Z"_M4X 

CI START = ^ . 

27r'f_Z"_MAX'R_START 

Wenn das Schalt kreisaquivalent mehrere RC-Glieder umfasst^ 
werden die fur die hoheren RC-Glieder erf orderlichen Start- 
werte fur die Fehlerminimierungsrechnung vorzugsweise aus 
den Dif f erenzspektren zwischen den gemessenen Daten und Da- 
ten ermittelt, die auf Grundlage der fur das erste RC-Glied 
berechneten Startwerte berechnet bzw. simuliert werden. 

Insbesondere die Auswahl sogenannter "guter" Startwerte, 
die nahe an den tatsachlichen Grofien der Bauelemente des 
Schalt kreisaquivalents liegen, verkurzen die Dauer der 
nach folgenden Fehlerminimierungsrechnung entscheidend . 
"Schlechte" Startwerte konnen hingegen dazu fuhren, dass 
die auf Grundlage der Startwerte durchgef uhrte Fehlermini- 
mierungsrechnung sinnlose Werte liefert- 

Bei einer bevorzugten Ausf iihrungsf orm des Verfahrens nach 
der Erfindung besteht das Schalt kreisaquivalent zur Simula- 
tion von Impedanzspektren auf Basis von bei einer Fehlermi- 
nimierungsrechnung gewonnenen Fitwerten aus einer seriellen 
Spaltung von vier RC-Gliedern. Durch Ermittlung von Start- 
werten fiir die einzelnen Bauelemente wird dann anhand eines 
Schwellenwertes bestimmt, wie viele RC-Glieder bei den Si- 
mulat ionsrechnungen berucksichtigt werden. Der Schwellen- 
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wert ist vorzugsweise ein von einem Benutzer voreinstellba- 
rer Wert. Uberpruft wird das prozentuale Verhaltnis des Wi- 
derstands des ersten RC-Glieds zum Widerstand des aktuellen 
RC-Glieds n gemali der Formel 



RCn_START > Wert[%] • RC1_START. 

Die GroBe ."Wert" ist die von dem Benutzer veranderbare Gro- 
fte, die beispielsweise auf 10 % voreingestellt ist. Wenn 
das Argument nicht erfullt ist, werden die Startwerte fur 
die Bauteile des betreffenden RC~Glieds auf Werte gesetzt, 
die keinen Einfluss auf eine Simulation der Impedanzspekt- 
ren haben. Diese Werte werden bei der Fehlerminimierungs- 
rechnung konstant gehalten. 

Fernef wird bei dem Verfahren die Validierungsgrolie be- 
stimmt, die die Ubereinstimmung zwischen dem errechneten, 
theoretischen Impedanzspektrum und dem jeweils zugeordne- 
ten, gemessenen Impedanzspektrum bewertet. 



Die Auswertegrolie stellt den fur die jeweilige Analyse re- 
levanten Ausgabewert dar, der beispielsweise bei der Er- 
mittlung eines Sensormaterials fur einen Gassensor eine 
Sensitivitat der jeweiligen Materialprobe wiedergibt. Die 
Sensitivitat ist ein MaS> fur die Gute eines Sensors. 

Die Sensitivitat eines resistiven Gassensors kann auf ver- 
schiedene Art und Weise definiert sein. Berucksichtigt man 
beispielsweise die Richtung der Anderung der Widerstande 
bei einer Gasbeauf schlagung, kann eine Sensitivitat S als 
Quotient des Widerstands R_TEST unter einer Prufgasatmo- 
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sphare und des Widerstands R_0 unter Ref erenzbedingungen 
wie folgt ausgedruckt werden 



^ R TEST ^ ^ 

^= =z J f^j2 oxidierende Gase 

R 0 



R 0 

S=z-{ = : fiir reduzierende Gase 

R TEST 



Das Vorzeichen der Sensitivitat S liefert eine Information 
hinsichtlich der Widerstandsanderung , 

Alternativ kann die Sensitivitat als Anderung der Wider- 
stande beschrieben werden, und zwar gemaB den Forineln 



^ , R TESTER 0 ^ 

S A = = fur oxidierende Gase 

R__TEST 

^ , R 0-R TEST ^ , . , 

o A = H = = ; fur reduzierende Gase- 

7? 0 



Bei diesen Formeln ergeben sich Sensitivitaten S_A zwischen 
-1 und 1, d- h. normierte Sensitivitaten. 



Es ist auch moglich, die Sensitivitaten beider Definitionen 
ineinander zu iiberfuhren. Die iiber den Wert S ausgedruckte 
Sensitivitat ist besonders aussagekraf tig bei einer groJien 
Anderung des Widerstands einer Materialprobe infolge einer 
Beauf schlagung mit einem Prufgas. Die iiber S_A ausgedruckte 
Sensitivitat ist besonders aussagekraf tig, wenn sich der 
Widerstand einer Materialprobe infolge einer Pruf gasbeauf - 
schlagung nur gering andert . Die durch S_A ausgedruckte 
Sensitivitat hat jedoch den Vorteil einer gro/ien Toleranz 
gegenuber Messungenauigkeiten, einer prazisen Beleuchtung 
eines Bereichs kleiner Widerstandsanderungen und damit ei- 
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ner besseren Bewertung von Querempf indlichkeiten sowie der 
Moglichkeit einer automatisierten Visualisierung und Daten- 
verarbeitung . 

Bei einer besonders einfachen Ausf uhrungsf orm des Verfah- 
rens nach der Erfindung besteht das Schalt kreisaquivalent 
aus einer virtuellen Anordnung realer elektronischer Bau- 
teile wie Kondensatoren und Widerstanden . In seiner ein~ 
fachsten Ausfiihrung besteht das Schalt kreisaquivalent aus 
einem Widerstand. 

Durch die Fehlerminimierungsrechnung, bei der eine Anglei- 
chung der Groiien der Bauelemente an die Messdaten erfolgt, 
wird ein simuliertes Impedanzspektrum fur das Schaltkreis- 
aquivalent fur die jeweilige Materialprobe ermittelt. Aus 
den angepassten GroBen der Bauelemente kann auf die elekt- 
rischen Eigenschaf ten der durch das Schaltkreisaquivalent 
beschriebenen Prozesse in der Materialprobe geschlossen 
werden. Wenn das elektrische Verhalten der Materialprobe 
durch mehrere Prozesse unterschiedlicher Relaxationszei ten 
bestimmt wird;. ist es erf orderlich, kompliziertere Schalt- 
kreisaquivalente heranzuziehen, die beispielsweise aus meh- 
reren in Reihe geschalteten RC-Gliedern bestehen. Wenn die 
Prozesse unterschiedlich auf eine Variation der Messbedin- 
gungen reagieren, kann eine Zuordnung einzelner Prozesse zu 
Bauelementen bzw. Bauelementgruppen vorgenommen werden, um 
so die einzelnen Prozesse separiert analysieren zu konnen. 
Durch die Fehlerminimierungsrechnung reduziert sich das Da- 
tenvolumen zur Beschreibung der Impedanzmessung auf die An- 
ordnung und die Grolien der Bauelemente. 
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Bei einer Probenplatte mit beispielsweise 64 Materialproben 
und einer Messung unter elf verschiedenen Gasatmospharen 
bei jeweils vier verschiedenen Temperaturen liegt das Mess- 
datenauf kommen bei 2816 Impedanzspektren . Die einzelnen Im- 
pedanzspektren werden gema/5 dem Verfahren nach der Erfin- 
dung simuliert und jeweils auf die abgeleiteten MessgroJien, 
die durch die Grofien der Bauelemente des Schalt kreisaquiva- 
lents wiedergegeben werden, reduziert. 

Die Berechnung eines theoretischen Impedanzspektrums eines 
Schaltkreisaquivalents erfolgt bei einer bevorzugten Aus- 
fuhrungsform des Verfahrens nach der Erfindung derart, dass 
eine komplexe Admittanz Y'^*^ und eine Phasenver schiebung cp 
der einzelnen Bauelemente der RC-Glieder bei einer gegebe- 
nen Winkelf requenz q (2n • Messf requenz ) gemaJi folgenden For- 
meln bestimmt wird: 



Durch eine Transformation konnen die jeweiligen Impedanzen 
ermittelt werden: 




(p = arctan(tyi?C) 



y = cos^ 



7" 



1 1 



r-+}' 



wobei Z' der Realteil der Impedanz und Z" der Imaginarteil 
der Impedanz ist. 
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Bei einer seriellen Schaltung der RC-Glieder konnen die Im- 
pedanzen der einzelnen RC-Glieder direkt summiert werden. 
Eine Berechnung der Impedanzen fur Frequenzen, die den 
Messf requenzen entsprechen, ergibt einen mit den Messdaten 
korrespondierenden Datensatz, so dass eine Fehlerabschat- 
zung zwischen dem gemessenen Datensatz und dem theoretisch 
ermittelten Datensatz moglich ist. 

Die bei dem Verfahren nach der Erfindung durchgef iihrte Feh- 
lerminimierungsrechnung wird bei einer vorteilhaf ten Aus- 
fuhrung durch Variation der Grofien der einzelnen Bauelemen- 
te um 1 % durchgef iihrt . Durch Analyse der Differenzen des 
theoretisch errechneten und des gemessenen Spektrums kann 
ein Fehler bestimmt werden. Wenn sich der Fehler nach einer 
Variation verkleinert, wird eine erneute Variation fur das 
gleiche Bauelement des Schalt kreisaquivalents durchgefuhrt 
und auf Grundlage dieser Variation ein theoretisches Impe- 
danzspektrum berechnet, mit dem eine neue Fehlerrechnung 
durchgefuhrt wird. Wenn sich der Fehler nicht verkleinert, 
wird die erfolgte Variation rlickgangig gemacht und die Va- 
riation der Bauelemente mit umgekehrten Vorzeichen variiert 
Oder ein anderes Bauelement variiert. 

Die Fehlerrechnung wird vorzugsweise so durchgefuhrt, dass 
die Funktion zur Bestimmung des Fehlers den Startwert des 
Widerstands des ersten RC-Glieds berilcksicht igt . Wenn der 
Wert groBer als ein nomineller Messwiders tand des bei der 
Me s sung eingeset zten Impedanz-Analysators , beispielsweise 
groBer als 3-10^ Q ist, wird der Fehler nur aus dem Imagi- 
narteil der berechneten Impedanz bestimmt. Um den hochfre- 
quenten Bereich des ermittelten Spektrums wichten zu kon- 
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neri;. warden die ermi ttelten Abweichungen vorzugsweise mit 

dem Logarithmus der Messfrequenz mult ipliziert . Diese Wich- 

tung ermoglicht es, physikalisch sinnlose oder auch fehler- 

hafte Messdaten im niederf requenten Bereich zu unterdru- 
cken . 



Zur Berechnung der Validierungsgroibe, d. h. zur Validierung 
der Gute des fur die jeweilige Materialprobe theoretisch 
berechneten Impedanzspektrums wird vorzugsweise ein Korri- 
dor um das theoretische Impedanzspektrum ermittelt, der 
beispielsweise 90 % der gemessenen Daten beinhaltet. Zur 
Minimierung der zur Berechnung der ValidierungsgroBe erfor- 
derlichen Zeit kann hierzu ein sukzessiver Approximations- 
Algorithmus genutzt werden, wobei die Grenzwerte des Algo- 
rithmus 0 Q und der zweifache Werte des Realteils der sum- 
mierten Impedanz bei der kleinsten Messfrequenz sein kon- 
nen . 



Das erf indungsgemafte Verfahren wird vorzugsweise zur Analy- 
se von Materialproben unter verschiedenen Pruf gasatmospha- 
ren und insbesondere bei verschiedenen Temperaturen einge- 
setzt. In diesem Fall werden fiir alle auf der Probenplatte 
angeordnete Materialproben Impedanzspektren unter den ver- 
schiedenen Messbedingungen gemessen . 

Die jeweils gemessenen Spektren werden dann bei dem Verfah- 
ren nach der Erfindung jeweils durch eine Fehlerminimie- 
rungsrechnung unter Zugrundelegung eines Schalt kreisaquiva- 
lents auf theoret ischem Wege simuliert. Dies resultiert in 
einer grolien Menge an Auswertegrofien, die die Zielgrolien 
bei dem Verfahren darstellen. 
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Die Zielgrof^en bzw. AuswertegroJien werden bei dem Verfahren 
nach der Erfindung vorzugsweise in eine Datenbank geschrie- 
ben und mittels einer Datamining-Funktionalitat ausgewer- 
tet- 

Bei der Datenbank, die zur Aufnahme und Bereitstellung von 
Datensatzen dient, handelt es sich zweckmafiigerweise um ei- 
ne relationale Datenbank, in der Inf ormationen nach Themen- 
kreisen geordnet in Form von Tabellen abgelegt sind. Bezie- 
hungen zwischen den einzelnen Datensatzen in den Tabellen 
werden durch sogenannte Identif ikationsschliissel herge- 
stellt - 

Die Datenbank kann beispielsweise weitere Eigenschaf ten der 
Materialproben, wie deren Synthesebedingungen der Edukte, 
deren Probengeschichte und dergleichen enthalten. Diese Ei- 
genschaften sind uber Tabellenbeziehungen miteinander ver- 
knupf t - 

Das Datamining kann mittels einer gegebenenf alls mehrdimen- 
sionalen Zielfunktion und/oder mittels einer visuellen Da- 
tamining-Funktionalitat durchgefuhrt werden. Das mittels 
einer Zielfunktion durchgef uhrte Datamining ist ein numeri- 
sches Verfahren, das auf den einzelnen Auswertegrolien ba- 
siert, die in der Datenbank beispielsweise auf geschlusselt 
nach Messtemperatur und Priifgas abgelegt sind. Dabei wird 
zunachst definiert, welche Eigenschaf ten von dem gesuchten 
Material verlangt werden. Beispielsweise kann bei der Suche 
nach einem Sensormaterial fur einen Gassensor angegeben 
werden, hinsichtlich welchen Gases das Sensormaterial sen- 
sitiv sein soil und welche Querempf indlichkeiten storen 
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konnten. Daraus ergibt sich also ein Anf orderungsprof il an 
einen Fingerabdruck der Sensitivitaten . 



Die visuelle Datamining-Funktionalitat arbeitet vorteilhaft 
derart, dass die Auswertegrolien der Materialproben der Pro- 
benplatte auf geschlusselt beispielsweise nach Prufgasen und 
Temperaturen dargestellt sind. 

Bei Einsatz der Vorrichtung nach der Erfindung und des Ver- 
fahrens nach der Erfindung liegt ein vorzugsweise vollauto- 
matisches High-Throughput-Impedanzsystem vor, mit dem neue 
Materialien mit einem hohen Probendurchsat z und einem ge- 
ringen Zeit- und Kostenauf wand entwickelt werden konnen. So 
konnen beispielsweise durch Einsatz des Systems an zwei Ta- 
gen 64 unterschiedliche Materialproben bei vier verschiede- 
nen Temperaturen und unter elf verschiedenen Prufgasatmo- 
spharen hinsichtlich ihrer sensorischen Eigenschaf ten un- 
tersucht und auch ausgewertet werden. 



Insbesondere kann das System zur allgemeinen Materialent- 
wicklung und speziell fur die Entwicklung von Sensoren im 
Kraftf ahrzeugbereich und im Sicherheitstechnikbereich ein- 
gesetzt werden. 



Die Erfindung hat auch eine Datenverarbeitungsanlage mit 
einem Datenverarbeitungsprogramm zur Durchfuhrung des er- 
f indungsgemafien Verfahrens zum Gegenstand. 

Weitere Vorteile und vorteilhafte Ausgestaltungen des Ge- 
genstandes nach der Erfindung ergeben sich aus der Be- 
schreibung, der Zeichnung und den Patentanspruchen . 
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Zeichnung 



Ein Ausf uhrungsbeispiel der Vorrichtung nach der Erfindung 
ist in der Zeichnung schematisch vereinfacht dargestellt 
und wird in der nachf olgenden Beschreibung in Verbindung 
mit einem Verfahren nach der Erfindung detailliert erlau- 
tert. Es zeigen 

Figur 1 einen Messaufbau einer Vorrichtung nach der Erfin- 
dung; 

Figur 2 eine Messeinrichtung der Vorrichtung nach Figur 1; 
Figur 3 einen Messkopf der Messeinrichtung nach Figur 2; 
Figur 4 eine Gas verteilungseinrichtung des Messkopfes nach 
Figur 3; 

Figur 5 einen Messablauf anhand eines Flussdiagramms ; 

Figur 6 ein Diagramm, das einen zeitlichen Messablauf fur 
Materialproben bei verschiedenen Temperaturen und unter 
verschiedenen Gasatmospharen darstellt; 

Figur 7 ein Diagramm, in dem ein Einlauf verhalten von vier 
unterschiedlich oberf lachendotierten Materialproben darge- 
stellt ist; 

Figur 8 ein Diagramm, in dem ein Ansprechzeitverhalten von 
drei verschiedenen Materialproben wahrend eines Prufgaspul- 
ses dargestellt ist; 

Figur 9 ein Schalt kreisaquivalent ; 

Figur 10 ein Flussdiagramm einer Fehlerminimierungsrech- 
nung; 

Figur 11 einen beispielhaf t en Verlauf einer Sensitivitat 
einer Materialprobe ; 

Figur 12 einen exemplar i schen Fingerabdruck von Sensitivi- 
taten einer Materialprobe hinsichtlich verschiedener Gasat- 
mospharen; und 
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Figur 13 eine stark schemat isierte Darstellung einer visu- 
ellen Datamining-Funktionalitat . 

Beschreibung des Ausfuhrungsbei spiels 

In den Figuren 1 bis 4 ist ein Messaufbau einer Vorrichtung 
10 zur Analyse einer Probenplatte 12 dargestellt, auf der 
64 Materialproben 13 angeordnet sind. Die Vorrichtung 10 
umfasst einen Hochtemperaturreaktor 14 r eine Gasversor- 
gungseinheit 16 und eine Mess- und Auswerteeinheit 18. 

Der Hochtemperaturreaktor 14, der insbesondere in den Figu- 
ren 2 bis 4 naher dargestellt ist, umfasst ein Gestell bzw. 
Gehause 20, in welchem ein hohenverstellbar.er Heizblock 22 
angeordnet ist, der an einer Ftihrungsstange 23 und an einer 
Gewindestange 24 gelagert ist, welche mit einer Kurbel 25 
versehen ist. 

Der Heizblock 22, der beispielsweise aus vier einen Heiz- 
raum begrenzenden Heizplatten mit einer Heizleistung von 
jeweils 1100 W besteht, ist so hohenverstellbar , dass der 
Plattentrager 28 und damit die Probenplatte 12 in den Heiz- 
raum eintauchbar ist. 

Des Weiteren umfasst der Hochtemperaturreaktor 14 einen 
Messkopf 26, der in einen mit dem Gestell 20 verbundenen 
Messkopf trager 27 eingesetzt ist und zur Kontakt ierung der 
64 auf der Probenplatte 12 angeordneten Materialproben 13 
mit der Mess- und Auswerteeinheit 18 dient. Die Probenplat- 
te 12 ist auf einem mit dem Messkopf trager 27 verbundenen 
Plattentrager 28 angeordnet. 
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Der Messkopf 26, der insbesondere in Figur 3 dargestellt 
ist, umfasst eine Grundplatte 29, die beispielsweise aus 
einer bearbeitbaren Glaskeramik gebildet ist und zur Halte- 
rung von 128 Messdrahten 30A, SOB dient, welche jeweils aus 
Platin bestehen und von einem Aluminiumoxidrohr uiimantelt 
sind. An ihren unteren Enden haben die Messdrahte 30A, 30B 
jeweils eine Schmelzkugel 31A, 31B, die zur Kontaktierung 
des jeweiligen Messdrahts an einer Kontaktf lache der Pro- 
benplatte 12 dient. Die Probenplatte 12 hat pro Material- 
probe zwei Kontaktf lachen, d. h. im vorliegendsten Fall 
insgesamt 128 Kontaktf lachen, die jeweils mit einem der 
Messdrahte 30A, 30B zusammenwirken . Die Messdrahte 30A, 30B 
sind an ihren oberen Enden im Bereich der Grundplatte 2 9 
jeweils mit einem vergoldeten Federkontakt 32A, 32B verse- 
hen, der bei fixiertem Messkopf 26 einen konstanten An- 
pressdruck der Messdrahte 30A, 30B bzw. der Schmelzkugel 
31A, 31B auf den jeweils zugeordneten Kontaktf lachen der 
Probenplatte 12 gewahrleistet . Die Ummantellungen der aus 
Platin gefertigten Messdrahte 30A, SOB sind des Weiteren an 
Halteplatten S3 und 34 fixiert, die an einem Gestange 35 
befestigt und parallel zu der Grundplatte 29 ausgerichtet 
sind . 



Die Federkontakte 32A, 32B sind iiber Leitungen 36A, 36B mit 
in dem Messkopf 26 integrierten SMB-Buchsen 37A, 37B ver- 
bunden, die wiederum uber abgeschirmte SMB-Leitungen 38A 
bzw. S8B mit der Mess- und Auswerteeinheit 18 der Vorrich- 
tung 10 verbunden sind. Insgesamt 128 SMB-Buchsen 37A, 37B 
hat der Messkopf 26, von denen jede mit einem Messdraht 
30A, 30B verbunden ist und an die jeweils eine zu der Mess- 
und Auswerteeinheit 18 fuhrende SMB-Leitung 38A, 388 ange- 
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schlossen ist. Der Ubersicht lichkeit halber sind in den Fi- 
guren 1 bis 3 aber jeweils nur zwei SMB-Buchsen 31A, 37B 
und zwei SMB-Leitungen 38A, 38B sowie die jeweils zugeord- 
neten Messdrahte 30A, 30B dargestellt. 

Der Messkopf 26 hat des Weiteren eine beispielsweise aus 
Quarzglas gefertigte, im Wesentlichen glockenf ormige Gas- 
verteilungseinrichtung 39, die in Figur 4 detailliert dar- 
gestellt ist und die uber eine Gasversorgungsleitung 40 aus 
Edelstahl mit der Gasversorgungseinheit 18 verbunden ist. 
Die glockenf ormige Gasverteilungseinrichtung 39 begrenzt 
einen Gasraum, der oberhalb der auf der Probenplatte 12 an- 
geordneten 64 Materialproben 13 angeordnet ist. 

Die Materialproben 13 sind beispielsweise aus Zinnoxid Sn02 
gebildet und weisen unterschiedliche Dotierungen auf, die 
beispielsweise aus Lanthaniden gebildet sind. Die Material- 
proben 13 sind matrixartig in acht Reihen und achte Spalten 
auf der Probenplatte 12 verteilt. 

Die in Figur 4 detailliert dargestellte Gasverteilungsein- 
richtung 39 hat in ihrem Gasraum 41 des Weiteren einen Dif- 
fusoreinsatz 42, der aus einer Quarzkugel gebildet ist und 
eine Vielzahl von Bohrungen 43 aufweist, die einen Durch- 
messer von jeweils etwa 1 mm haben . 

Urn einen gleichmaliigen Gasaustritt aus dem Gasraum 41 zu 
gewahrleisten, hat die glockenf ormige Gasverteilungsein- 
richtung 39 an ihren Randern Distanzhalter 43, die einen 
Spalt von 0,8 mm Breite zwischen der Gasverteilungseinrich- 
tung 39 und der Probenplatte 12 festlegen. 
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Zur Beauf schlagung des Gasraums 41 und damit der Material- 
proben 13 mit unter schiedlichen Priif gasatmospharen weist 
die Gasversorgungseinheit 16 zwei hier nicht naher darge- 
stellte Gasf laschenschranke mit je vier Gasflaschen auf, 
die jeweils mit einem Gasdurchf lussregler 44B, 44C, 

44D, 44E, 44F, 44G bzw. 44H verbunden sind, wobei in einer 
Gasflasche feuchte synthetische Luft;. in der zweiten Gas- 
flasche Wasserstoff, in der dritten Gasfalsche Methan, in 
der vierten Gasflache synthetische Luft, in der fiinften 
Gasflache Stickstof f dioxid, in der sechsten Gasflasche 
Stickstof fmonoxid, in der siebten Gasflache Propen und in 
der achten Gasflaschen Kohlenstof f monoxid enthalten ist. 
Die Kapazitaten der Durchf lussregler 44A, 44B, 44C, 44D, 
44E und 44F liegen jeweils zwischen 0 seem und 100 seem. 
Die Kapazitaten der Durchf lussregler 44G und 44H liegen je- 
weils zwischen 0 seem und 10 seem. Uber entsprechende Steu- 
erung des Volumenstroms der verschiedenen Case konnen mit- 
tels der Batterie aus den acht Gasdurchf lussreglern 44A, 
44B, 44C, 44D, 44E, 44F, 44G und 44H Prufgase unterschied- 
licher Zusammenset zungen in eine Sammelleitung 45 einge- 
speist warden, die mit der Gaszuf uhrleitung 40 verbunden 
ist . 

Urn eine relative Feuchtigkeit des jeweiligen Prufgases ein- 
zustellen, kann dem Priifgas ein feuchtes Tragergas, das 
beispielsweise aus synthet ischer Luft besteht, zugemischt 
warden. Die Feuchte des Tragergas wird dadurch eingestellt, 
dass es durch ein Wasserreservoir 46 geleitet wird. Eine 
Messung des angef euchtet en Tragergases mittels eines hier 
nicht naher dargestellten Feuehtesensors ergibt beispiels- 
weise eine relative Feuchte von etwa 90 % bei Raumtempera- 
tur . 
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Die Mess- und Auswerteeinheit 18 umfasst zwei Relaisschalt- 
felder 50 und 51, an die jeweils 64 Leitungen 38A bzw. 38B 
angeschlossen sind, die zu dem Messkopf 26 des Hochtempera- 
turreaktors 14 fuhren. Die Relaisschaltf elder 50 und 51 
bilden jeweils eine 3x64-Matrix aus hochf requenztauglichen 
Relais . 

Die Relaisschaltf elder 50 und 51 sind uber eine digitale 
Steuerleitung 52 mit einem Mess- und Auswerterechner 53 
verbunden. Des Weiteren sind die Relaisschalt f elder 50 und 
51 liber Messleitungen 54 mit einem Impedanzanalysator 64 
und einem sogenannten Source-Meter 55 verbunden. Diese bei- 
den Messgerate sind ebenfalls uber die digitale Steuerlei- 
tung 52 mit dem Mess- und Auswerterechner 53 verbunden. Die 
Adressierung des Impedanzanalysators 54 und des Source- 
Meters 55 erfolgt iiber die beiden Relais-Schaltf elder 50 
und 51. 

Der Mess- und Auswerterechner 53 ist ferner uber eine wei- 
tere digitale Steuerleitung 56 mit einem D/A-A/D-Wandler 57 
verbunden, der uber eine analoge Steuerleitung 58 zum einen 
mit der Heizeinrichtung 22 des Hochtemperaturreaktors 14 
und zum anderen mit den Gasdurchf lussreglern 44A, 44B, 44C, 
44D, 44E, 44F, 44G und 44H der Gasversorgungseinheit bzw. 
Gasmischbatterie 16 verbunden ist. 

Auf dem Mess- und Auswerterechner 53 ist eine modulare 
Mess- und Steuersof tware abgelegt, die uber eine Skript- 
steuerung eine vollstandige Automat isierung von mittels der 
Vorrichtung 10 durchgef uhrter Messungen ermoglicht. Durch 
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den modularen Aufbau der Mess- und Steuersof tware ist eine 
problemlose Erweiterung des Messsystems moglich. 

Mittels des Source-Meters 55 konnen Gleichstromwiderstande, 
U/I-Kennlinien oder Spannungen der einzelnen Materialproben 
13 gemessen werden. Diese Werte sowie die mittels des Impe- 
danz-Analysators 54 ermittelten Messwerte konnen uber die 
digitale Steuerleitung 52 an den Mess- und Auswerterechner 
53 ubergeben werden, welcher eine Datenbank fur die Messda- 
ten umfasst. 

Zur Vermessung der auf der Probenplatte 12 angeordneten Ma- 
terialproben 13 in einem Hochdurchsat z-Modus bzw. nach ei- 
nem High-Throughput-Verf ahren wird eine Skript sprache ein- 
gesetzt, wobei uber eine Skriptdatei eine Liste von Aufga- 
ben, die aus Schlusselwortern und Parametern bestehen, an 
die Software ubergeben wird, mittels welcher die Skriptda- 
tei abgearbeitet wird und welche samtliche Funktionen des 
Systems steuert. Die Skriptdateien gewahrleisten eine kon- 
tinuierliche Uberprufung von Regelparametern, so dass Mes- 
sungen unter falschen Messbedingungen dadurch ausgeschlos- 
sen werden, dass eine weitere Verarbeitung der Skriptdatei 
ausgesetzt wird. Die Skriptdatei ist hinsichtlich ihrer 
Lange nicht begrenzt. 



Ein mittels der Vorrichtung nach Figur 1 durchf uhrbarer 
Messablauf ist in Figur 5 anhand eines Flussdiagramms dar- 
gestellt. Zur Analyse von resistivgassensorischen Eigen- 
schaften werden die Materialproben 13 der Probenplatte 12 
zum einen unter einer Ref erenzgasatmosphare, die beispiels- 
weise aus synthetischer Luft einer relativen Feuchtigkeit 
von beispielsweise 45 % gebildet ist, und unter verschiede- 
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nen Pruf gasatmospharen elektrisch charakter isiert . Durch 
eine Modulation der Messtemperatur konnen Inf ormationen 
uber den Einfluss der jeweiligen Operations temperatur auf 
die als Sensormaterial geeigneten Materialproben sowie uber 
die Aktivierungsenergie von Leitf ahigkeit sprozessen gewon- 
nen werden. 



Urn eine Vergleichbarkeit von Messdaten mehrerer Probenplat- 
ten zu gewahrleisten, ist die eingesetzte Skriptdatei als 
Standardskript aufgebaut, welches das komplette Hochdurch- 
satz-Screening steuert, dessen chronologischer Ablauf in 
dem in Figur 5 dargestellten Flussdiagramm beschrieben ist. 
Hierbei wird in einem ersten Ver f ahrensschr itt Ml eine 
Messtemperatur T_MESS der Probenplatte 12 eingestellt. Mes- 
sungen konnen beispielsweise in einem Temperaturbereich 
zwischen 400 °C und 250 °C in Schritten von 50° durchge- 
fuhrt werden. Wahrend der Abkuhlphasen wird die Probenplat- 
te 12 mit einem Referenzgas aus synthetischer Luft einer 
relativen Feuchtigkeit von 45 % beauf schlagt . 

Bei einer Anderung der Probentemperatur stellt sich ein 
stabiler Grundwiderstand der Materialproben erst nach einer 
bestimmten Zeit ein, so dass eine Konditionierung der Mate- 
rialproben in einem Schritt M2 erforderlich ist. Tempera- 
turanderungen der Materialproben fuhren zu einem metastabi- 
len Zustand intrinsischer Defekte, die Sauerstof f leerstel- 
len darstellen konnen und deren thermodynamische Gleichge- 
wichtseinstellung endliche Zeit benotigt. In Figur 7 ist 
beispielhaft das Konditionierungs- bzw. Einlauf verhalten 
fur vier unterschiedlich oberf lachendot ierte In203-Proben 
bei Erreichen einer Zieltemperatur von 300 °C dargestellt. 
Die Widerstande R der Proben nahern sich mit der Zeit 
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asymptotisch einem Grenzwert, der den sogenannten Grund- 
bzw, Ref erenzwiderstand darstellt. Dieser Grenzwert wird 
extrapoliert nach etwa 90 Minuten erreicht. Zur Gewahrleis- 
tung eines konstanten Ref erenzwiderstands 1st es vorteil- 
haft, eine Konditionierungszeit von 120 Minuten zu wahlen. 



AnschlieJiend wird das zur Messung erf orderliche Prufgas in 
einem Verf ahrensschritt M3 zusammengestellt und in den uber 
der Probenplatte 12 angeordneten Gasraum 41 eingeleitet, 
Das erste Prufgas umfasst beispielsweise Wasserstoff mit 
einer Konzentration von 25 ppm. Als Tragergas dient synthe- 
tische Luft. Der Gasstrom wird mit feuchter synthet ischer 
Luft auf 100 seem eingestellt. 

Urn zu gewahrleisten, dass die einzelnen Materialproben un- 
abhangig von ihrer relativen Position auf der Probenplatte 
ihren Grundwider stand erreicht haben, erfolgt anschlieJiend 
in einem Schritt iyi4 eine Vorlaufphase fur das Prufgas, die 
zuvor untersucht wurde. In Figur 8 ist der Widerstandsver- 
lauf von drei auf einer Diagonalen der Probenplatte liegen- 
den Materialproben wahrend eines 40 Minuten andauernden 
Pruf gaspulses dargestellt, wobei das Pulsgas Propen einer 
Konzentration von 50 ppm umfasst. Das Basismaterial der Ma- 
terialproben ist bei dem in Figur 8 dargestellten Beispiel 
wegen seiner hohen Sensitivitat auf Kohlenwasserstof f e 
Zinnoxid Sn02 . Unabhangig von der relativen Position auf 
der Probenplatte sinkt der Widerstand der Materialproben 
aufgrund des Pruf gaspulses innerhalb von etwa 6 Minuten auf 
einen konstanten Wert. Die Ansprechzeit betragt also je- 
weils 6 Minuten, wobei das Ansprechverhalten im Wesentli- 
chen unabhangig von der Position der jeweiligen Material- 
probe auf der Probenplatte ist. Nach Beendigung des Pruf- 



wo 2005/061093 



27 



PCT/DE2004/002539 



gaspulses nahern sich die Widerstande innerhalb von etwa 10 
Minuten asymptotisch den Grundwiderstanden an. Das An- 
sprechverhalten ist demnach im Wesentlichen unabhangig von 
der Position der Materialprobe auf der Probenplatte . Die 
fur die Skriptdatei gewahlte Vorlaufzeit jedes Priif- oder 
auch Ref erenzgases betragt vor einer entsprechenden Messung 
etwa 15 Minuten. 

In einem anschlieJienden Verf ahrensschritt M5 erfolgt die 
Messung von Impedanzspektren fur die 64 auf der Probenplat- 
te angeordneten Materialproben . Dabei werden in der Skript- 
datei fur die Messungen der Impedanzspektren als Parameter 
f es tgelegt : 



Amplitude der Messspannung [V] 0,1 

Startf requenz [Hz] 10 

Endfrequenz [Hz] 10^ 

Messpunkte pro Frequenzdekade 15 

Bias [V] 0 

Modus [HS: Highspeed; NO: Normal; AV: Average] HS 



Die gewonnenen Messdaten werden mittels des Impedanzanaly- 
sators 64 ermittelt und an den Mess- und Auswerterechner 53 
bzw. die darauf abgelegte Datenbank ubergeben. 

Im Anschluss daran kann die Messung unter einer anderen 
Priif gasatmosphare durchgefuhrt werden, wobei hierbei wie- 
derum auch die Schritte M3, M4 und M5 durchlaufen werden. 
Beispielsweise konnen als weitere Prufgase Kohlenstof f mono- 
xid einer Konzentration von 50 ppm, Stickstof f monoxid einer 
Konzentration von 5 ppm, Stickstof fdioxid einer Konzentra- 
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tion von 5 ppm oder auch Propen einer Konzentration von 
25 ppm eingesetzt werden. 



Das gesamte Screening kann dann durch Zuruckspringen zu dem 
Schritt Ml bei einer anderen Messtemperatur durchgefiihrt 
werden . 



In Figur 6 ist ein durch ein Standard-Skript bestimmter 
Messablauf visualisiert , wobei in diesem Diagramm Messungen 
unter Ref erenzbedingungen und wahrend Konditionierungspha- 
sen der Ubersichtlichkeit halber nicht dargestellt sind. 
Messungen unter Ref erenzbedingungen , d. h. unter einer Re- 
f erenzgasatmosphare, werden stets vor Einleiten eines neuen 
Prufgases H2, CO, NO, NO2 oder Propen bei der betreffenden 
Temperatur T durchgef uhrt . Es liegt stets jeweils eine so- 
genannte Priifgas-Vorlauf phase X sowie eine Messphase Y vor. 
Ferner sind in Figur 6 die jeweiligen Konzentrationen C der 
Pruf gase dargestellt . 



Die Mess- und Auswertesof tware kann eine Funktionalitat ha- 
ben, bei der das jeweils gemessene Impedanzspektrum gra- 
phisch an einem Monitor dargestellt wird. Die Messdaten, 
die vorzugsweise im ASCII-Format in der Datenbank oder in 
den den Probenplatten zugeordneten Verzeichnissen auf dem 
Mess- und Auswerterechner gespeichert werden, sind somit 
direkt einer visuelien Kontrolle zuganglich. 



Des Weiteren kann die Mess- und Auswertesof tware mit einer 
Funktionalitat ausgestattet sein, die Impedanzspekt ren der 
Probenmaterialien einer Probenplatte positionsabhangig in 
einer Bildmatrix darstellt. Dabei konnen sowohl Rohdaten 
Oder auch abgeleitete Daten dargestellt werden. Die Funkti- 
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onalitat kann auch ein Auswertef enst er umfassen^ in das Da- 
ten einer bestimmten Materialprobe iibertragen werden kon- 
nen. Auch kann die Funktionalitat mit weiteren bildgebenden 
Mess- Oder Auswertesystemen kombiniert werden. Bilddaten 
von diesen Systemen konnen dann ebenfalls positionsabhangig 
in die Bildmatrix eingelesen werden, urn so weitere Informa- 
tionen uber die Proben zu erhalten. 



Nach Messung der Impedanzspektren fur die einzelnen Materi- 
alproben wird fiir jede Materialprobe ein theoretisches Im- 
pedanzspektrum auf Basis eines Schalt kreisaquivalent s be- 
rechnet . Das Schalt kreisaquivalent besteht im vorliegenden 
Fall aus einer seriellen Schaltung von vier RC-Gliedern der 
in Figur 9 dargestellten Art. Die theoretisch berechneten 
Impedanzspektren werden dabei jeweils bestmoglich an das 
korrespondierend gemessene Impedanzspektrum angepasst, und 
zwar dadurch, dass die Grofien bzw. Dimensionierungen der 
einzelnen Bauelemente des Schaltkreisaquivalent s variiert 
werden . 



Die durchgef uhrte Fehlerminimierungsrechnung wird nachfol- 
gend anhand des in Figur 10 dargestellten Flussdiagramms 
erlautert . 



Zur Durchfiihrung der Fehlerminimierungsrechnung ist es er- 
forderlich, fiir die einzelnen Bauelemente des Schaltkreis- 
aquivalents Startwerte zu bestimmen. Gleichzeitig ist es 
erf orderlich, die Anzahl der bei der Fehlerminimierungs- 
rechnung zu berucksicht igenden RC-Glieder zu bestimmen. Zu 
diesen Zwecken wird in einem Schritt SI die maximale imagi- 
nare Impedanz Z"_MAX aus den Messdaten fiir die betreffende 
Materialprobe ermittelt . Bei der Auswahl wird auch beur- 
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teilt, ob es sich bei der maximal gemessenen imaginaren Im- 
pedanz um eine Fehlmessung handelt, was durch eine Untersu- 
chung der bei benachbarten Messpunkten gemessenen imagina- 
ren Impedanzen, d. h. durch eine Untersuchung des lokalen 
Maximums , er f olgt . 

Anschlieftend warden in einem Schritt S2 Startwerte R1_START 
fur den Widerstand und C1_START fur die Kapazitat des ers- 
ten RC-Glieds auf der Basis von Z"_MAX und der korrespon- 
dierenden Messf requenz gemaB den Formeln 

R\_START^ -2-Z"_AdAX 

CI START = ? 

In • f _Z'' _A4AX ' R _START 

berechnet, wobei f_Z"_MAX die Messfrequenz bei der maximal 
gemessenen, imaginaren Impedanz ist. 

Anschlieftend werden in einer Schleife S3 Startwerte fur die 
Widerstande und Kapazitaten der weiteren drei RC-Glieder 
des Schaltkreisaquivalents ermittelt . 

Hierzu wird zunachst in einem Schritt S4 auf Basis der 
Startwerte R1_START und C1_START ein theoretisches Impe- 
danzspektrum berechnet und in einem Schritt S5 ein Diffe- 
renzspektrum zwischen dem theoretischen Impedanzspektrum 
und dem gemessenen Impedanzspektrum berechnet und aus die- 
sem Dif f erenzspektrum wiederum das Maximum der imaginaren 
Impedanz Z"_MAX ermittelt. Ausgehend von diesem Maximum 
werden dann in einem Schritt S6 die Startwerte Rn_START und 
Cn_START (n = 2 bis 4) fur das jeweils betrachtete RC~Glied 
gemaft den im Zusammenhang mit Schritt S2 erlauterten For- 
meln berechnet. 
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AnschlieiJend wird in einem Schritt S7 eine Schwellenwertbe- 
trachtung fur den berechneten Widerstandsstartwert Rn_START 
durchgef iihrt , mittels dem bestimmt wird, ob das jeweilige 
RC-Glied bei den anschlieiJenden Simulationsrechnungen be- 
rucksichtigt warden soil. Der Schwellenwert ist ein von dem 
Benutzer veranderbarer Wert. Uberpruft wird bei dem Ver- 
gleich das prozentuale Verhaltnis des Widerstands des ers- 
ten RC-Glieds zum Widerstand des aktuellen Bauteils. 

Wenn der berechnete Startwert Rn_START grolier ist als der 
Schwellenwert, wird das RC-Glied berucksichtigt , zu Schritt 
S3 zuruckgesprungen und eine Berechnung der Startwerte 
Rn__START und Cn_START fur das nachste RC-Glied vorgenommen . 
Wenn das Argument nicht erfullt ist, werden die Groften der 
Bauteile der hoheren RC-Glieder auf Werte festgelegt, die 
keinen Einfluss auf die Simulat ionsrechnung fur ein Impe- 
danzspektrum haben. Diese Werte bleiben bei nachf olgenden 
Anpassungsschritten konstant und betragen beispielsweise 
Rn_START=l und Cn_START=10'^\ Diese Festlegung erfolgt in 
einem Schritt S8. 

Die Schwellenwertabf rage bestimmt also die Anzahl der RC- 
Glieder, die in den folgenden Schritten berucksichtigt wer- 
den. Wenn der Schwellenwert unterschritten wird, wird di- 
rekt zur Fehlerminimierungsrechnung auf der Basis der An- 
zahl der bestimmten RC-Glieder, beispielsweise auf der Ba- 
sis eines RC-Glieds ubergegangen, wobei hierbei eine 
Schleife S9 uber alle m berucksichtigten RC-Glieder durch- 
gefuhrt wird. 

Innerhalb der Schleife S9 wird eine weitere Schleife SIO 
Lauf = 1 to 3 durchgef uhrt , wobei es sich hierbei urn eine 
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empirische Zahl von Durchlaufen handelt, die zur Erhohung 
der Genauigkeit der Fehlerminimierungsrechnung durchgefuhrt 
werden . 

Innerhalb der Schleife SIO wird zunachst in einem Schritt 

511 als Variable VAR der Widerstandsstartwert Rl START des 
ersten RC-Glieds herangezogen . Die Variation der einzelnen 
Bauelemente bei der Fehlerminimierungsrechnung betragt 1 %, 
was in einem Wert FAK = 0,01 in Schritt Sll zum Ausdruck 
kommt. Der Fehler wird auf 10^^ voreingestellt . Die Varia- 
tion des jeweiligen Bauelements erfolgt in einem Schritt 

512 gemali der Formel VAR = V-FAK. Fine Berechnung eines 
theoretischen Impedanzspektrums erfolgt dann auf Grundlage 
des variierten Werts in einem Schritt S13. 

AnschlieJiend wird ein Fehler berechnet, der auf einem* Ver- 
gleich zwischen dem fur die betreffend angepassten Bauele- 
mentgroBen berechneten Impedanzspektrum und dem gemessenen 
Impedanzspektrum beruht . Fine Funktion zur Bestimmung des 
Fehlers ist abhangig von dem Startwert des Widerstands des 
ersten RC-Glieds. Zur Berechnung des Fehlers wird in einem 
Schritt S14 uberprtift, ob der jeweilige Startwert groBer 
als der nominelle Messwiderstand des eingesetzten Impedanz- 
analysators, beispielsweise grofier als 3 x lO"^ Q ist. Ist 
dies der Fall, wird der Fehler in einem Schritt S15 nur 
durch Betrachtung der Imaginarteile der Impedanzen gemaB 
der Formel 

Fehler = |Z''_fit - Z''_mess| • logf 

bestimmt, wobei zur Wichtung des Hochf requenzbereichs des 
Spektrums eine Multiplikation mit dem Logarithmus der Mess- 
frequenz erfolgt, wodurch physikalisch sinnlose bzw. feh- 
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lerhafte Messdaten im Niederf requenzenbereich unterdruckt 
werden. Anderenfalls erfolgt in einem Schritt S16 eine 
Fehlerbetrachtung auf Grundlage sowohl der Realteile als 
auch der Imaginarteile der Impedanzen gemaft der Formal 

Fehler = |Z''_fit - Z ' '^mess | + | Z '__f it - Z'_mess|. 

In einem Schritt S17 wird uberpruft, ob der Fehler durch 
die Variation des Bauelements in Schritt S12 kleiner gewor- 
den ist. Wenn dies der Fall ist, wird das betreffende Bau- 
teil, beispielsweise der Widerstand des betreffenden RC- 
Glieds durch Zuruckspringen zu dem Schritt 312 erneut vari- 
iert- Wenn der Fehler nicht gesunken ist, wird die im 
Schritt 812 erfolgte Variation in einem Schritt S18 ruck- 
gangig gemacht und in einem Schritt S19 festgelegt, ob das 
Vorzeichen der Variation in einem Schritt S20 geandert oder 
ob in einem Schritt S21 und einem nachf olgenden Schritt S22 
das nachste Bauelement, in diesem Falle die Kapazitat C, 
als variables Bauelement gewahlt wird. In diesem Fall er- 
folgt dann abermals ein Rucksprung zu dem Schritt S12, wo- 
bei die Groiie des Bauteils dann wiederum so lange variiert 
wird, wie der Fehler sinkt. 

Die Variation des Widerstands und die Variation der Kapazi- 
tat wird durch Rucksprung zu dem Schritt SIO fur jedes RC- 
Glied zweimal wiederholt. Danach erfolgt in Abhangigkeit 
von der Anzahl der zu berucksichtigenden RC-Glieder ein 
n-maliger Rucksprung zu dem Schritt S9- 

Wenn keine weitere Minimierung des Fehlers durch Variation 
des Widerstands und der Kapazitat des RC-Glieds moglich 
ist, werden die den kleinsten Fehler ergebenden Werte fur 
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'den Widerstand und die Kapazitat als Auswertegroften der 
Fehlerminimierungsrechnung ausgegeben . 



Nach Abschluss der in den Schritten S9 bis S22 erfolgenden 
Fehlerminimierungsrechnungen wird fiir das theoretisch er- 
mittelte Impedanzspektrum in einem Schritt S23 ein Validie- 
rungsf unktion ermittelt, durch die die Gute des theoretisch 
berechneten Spektrums abgeschatzt werden kann. Bei der Be- 
stimmung der Validierungsf unktion wird zur Abschatzung ei- 
ner Gute der Fehlerminimierungsrechnung ein Korridor urn das 
errechnete Impedanzspektrum ermittelt, der 90 % der Messda- 
ten beinhaltet- Der so ermittelte Fehler liegt zwischen 0 
und 1, wobei ein Fehler von 0 eine ideale Ubereinstimmung 
zwischen der Messung und der Simulat ionsrechnung und ein 
Fehler von 1 keine Ubereinstimmung zwischen der Messung und 
der Simulationsrechnung zum Ausdruck bringt . 

Bei Einsatz des Verfahrens zur Festlegung eines Sensormate- 
rials eines Gassensors werden auf Grundlage der bei der 
Fehlerminimierungsrechnung gewonnenen Wider stande der 
Schaltkreisaquivalente Sensitivitaten S_A berechnet. Die so 
gewonnenen Sensitivitaten S_A sind normiert und liegen zwi- 
schen -1 und +1, wie dem Diagramm in Figur 11 zu entnehmen 
ist. In Figur 11 sind Sensitivitaten fur einen angenommenen 
Ref erenzwiderstand von 100 Q bei variierendem Testwider- 
stand R_test dargestellt. Die Sensitivitaten S_A werden in 
der Datenbank abgelegt. 



Anschlieftend wird auf Grundlage der ermittelten Sensitivi- 
taten S_A ein numerisches Datamining durchgef uhrt , um fiir 
einen speziellen Anwendungsf all einen optimalen Gassensor 
auszuwahlen. Das Datamining wird insbesondere auf numeri- 
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schem Wege durchgef iihrt , wobei auf die einzelnen Sensitivi- 
taten S_A zugegriffen wird, die in der Datenbank aufge- 
schlusselt nach Messtemperatur und Prufgas abgelegt sind. 
Zunachst wird definiert, welche Eigenschaf ten von dem ge- 
wunschten Sensormaterial verlangt werden. Im einfachsten 
Fall wird angegeben, fur welches Prufgas der Sensor einge- 
setzt werden soli und welche Querempf indlichkeiten storen. 
Daraus ergibt sich ein Anf orderungsprof il an den Fingerab- 
druck der Sensitivitaten, wie in Figur 12 dargestellt ist. 
Das Anforderungsprof il wird umgesetzt in sogenannte ">" 
bzw. "<" -Anforderungen, wodurch Datensatze mit den ge- 
wiinschten Eigenschaf ten charakterisiert werden. Weitere Be- 
dingungen, wie die Art des Basismaterials des Sensormateri- 
als Oder dessen Oberf lachendotierung konnen in die Ziel- 
funktion auf genommen werden. Bei dem Beispiel in Figur 12 
gestaltet sich das Anf orderungsprof il derart, dass alle 
Priifgase A, B, C, E aufter dem Prufgas D eine Sensitivitat 
kleiner 0,2 und groBer -0,2 haben und das Prufgas D eine 
Sensitivitat grower 0,9 hat. Mit X ist jeweils die Referenz 
angegeben. Das Anf orderungsprof il wird mittels einer Aus- 
wertefunktionalitat als SQL-Filterabf rage umgesetzt und di- 
rekt an die Datenbank iibergeben. Die Ergebnisse der SQL- 
Anweisung konnen dann in tabellarischer Form auf einem 
Bildschirm dargestellt werden. 

Alternativ oder zusatzlich konnen die fiir die verschiedenen 
Materialproben bei den verschiedenen Messbedingungen ermit- 
telten Sensitivitaten mittels einer visuellen Datamining- 
Funktionalitat ausgewertet werden. Bei dieser Funktionali- 
tat, die in Figur 13 dargestellt ist, werden Sensitivitaten 
sogenannter Bibliotheksplatten auf geschlusselt nach Prufga- 
sen und Temperaturen auf einem Bildschirm dargestellt. In 
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der beispielhaf ten Darstellung nach Figur 13 sind Biblio- 
theksplatten 101 bis 112 fur vier Temperaturen A, B, C und 
D und fur drei verschiedene Prufgase I, II und III darge- 
stellt. Jede Bibliotheksplatte 101 bis 112 ist einer Tempe- 
ratur A, B, C bzw. D und einem Prufgas I, II bzw. Ill zuge- 
ordnet. Die Sensitivitaten S_A der einzelnen Materialien 
der Bibliotheksplatten 101 bis 112 sind jeweils entspre- 
chend ihrer Position auf der Probenplatte als Kreise 120 
und in Fehlfarben dargestellt. In Figur 13 sind der Uber- 
sichtlichkeit halber fur jede Bibliotheksplatte 101 bis 112 
nur vier von 64 Materialproben dargestellt. Positive Sensi- 
tivitaten S_A werden beispielsweise in Farben dargestellt, 
die von Schwarz uber Rottone bis Gelb gehen, wohingegen ne- 
gative Sensitivitaten beispielsweise in Farben dargestellt 
werden, die von Schwarz uber Blautone bis Tiirkis gehen. Der 
Durchmesser der einzelnen Kreise 120 wird bestimmt durch 
die Validierung der Fehlerminimierungsrechnung, wobei bei 
einem groBen Fehler ein kleiner Kreis und bei einem kleinen 
Fehler ein vergleichsweise groBer Kreis dargestellt wird. 
Damit ist fiir einen Benutzer eine intuitive Beurteilung 
grofter Datenmengen moglich. 

Die zu der anhand Figur 13 dargestellten Funkt ionalitat er- 
forderlichen Datensatze werden direkt aus der Datenbank 
extrahiert, so dass uber Filter funkt ionen die darzustellen- 
den Messungen ausgewahlt werden konnen. Die Anzahl der an- 
gezeigten Datensatze ist nur durch den Arbeitsspeicher des 
eingesetzten Auswerterechners beschrankt- 

Des Weiteren ist es bei der visuellen Datamining-Funktiona- 
litat moglich, dass zusatzliche Informat ionen zu den ein- 
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zelnen Materialproben durch bildschirmgestut zte Auswahl ei- 
nes bestimmten Kreises angezeigt warden. 
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Anspriiche 

1. Vorrichtung zur Analyse einer Probenplatte (12), auf 
der mindestens zwei Materialproben (13) angeordnet 
sind, umfassend einen Trager*(28) fur die Probenplatte 
(12) und Kontaktiermittel zur elektrischen Kontaktie- 
rung der Materialproben (13), gekennzeichnet durch ei- 
nen in einem Gehausetrager (27) einsetzbaren Messkopf 
(26), der zur elektrischen Verbindung mit den Kontak- 
tiermitteln je Materialprobe (13) zwei Messdrahte (30A, 
30B) umfasst, die mit Vorspannung an Kontaktf lachen der 
Probenplatte (12) anliegen und mit einer Mess- und Aus- 
werteeinheit (18) verbunden sind. 

2- Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, 

dass die Messdrahte (30A, 30B) uber Schmelzkugeln (31A, 
31B) an den Kontaktf lachen der Probenplatte (12) anlie- 
gen. 

3. Vorrichtung nach Anspruch 1 oder 2, dadurch gekenn- 
zeichnet, dass die Messdrahte (30A, 30B) jeweils mit 
einem Federkontakt (32A, 328) verbunden sind, der einen 
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konstanten Anpressdruck des jeweiligen Messdrahts (30A, 
30B) auf der jeweiligen Kontaktf lache gewahrleistet , 

4. Vorrichtung nach einem der Anspriiche 1 bis 3, dadurch 
gekennzeichnet , dass der Messkopf (26) mit einer Gas- 
versorgungseinheit (16) verbunden ist. 

5. Vorrichtung nach Anspruch 4, dadurch gekennzeichnet , 
dass die Gasversorgungseinheit (16) mit einer Datenver- 
arbeitungseinheit (53) der Mess- und Auswerteeinheit 
(18) verbunden ist. 

6. Vorrichtung nach Anspruch 4 oder 5, dadurch gekenn- 
zeichnet, dass die Gasversorgungseinheit (16) eine Gas- 
mischeinrichtung umf asst . 

7. Vorrichtung nach einem der Anspruche 4 bis 6, dadurch 
gekennzeichnet , dass die Gasversorgungseinheit (16) ein 
Wasserreservoir (46) umfasst. 

8. Vorrichtung nach einem der Anspruche 4 bis 1, dadurch 
gekennzeichnet, dass der Messkopf (26) zur Gasbeauf- 
schlagung der Probenplatte (12) einen vorzugsweise von 
einer im Wesentlichen glockenf ormigen Verteilereinrich- 
tung (39) gebildeten Gasraum umfasst, der mit der Gas- 
versorgungseinheit (16) verbunden ist. 

9. Vorrichtung nach Anspruch 8, dadurch gekennzeichnet, 
dass in dem Gasraum ein Diffusor (42) angeordnet ist. 



10. 



Vorrichtung nach Anspruch 8 oder 9, dadurch gekenn- 
zeichnet, dass der Gasraum mit einem Gasauslass verse- 
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hen ist, der vorzugsweise durch mindestens einen Dis- 
tanzhalter (43) gebildet ist, der zwischen der Proben- 
platte (12) und der Verteilereinrichtung (39) angeord- 
net ist. 

11. Vorrichtung nach einem der Anspruche 1 bis 10, dadurch 
gekennzeichnet, dass die Mess- und Auswerteeinheit (18) 
zwei Relaisschaltfelder (50, 51) umfasst, die mit den 
Messdrahten (30A, 30B) verbunden sind und vorzugsweise 
jeweils eine 3x64-Matrix auch hochf requenztauglichen 
Relais aufweisen. 

12. Vorrichtung nach einem der Anspruche 1 bis 11, dadurch 
gekennzeichnet , dass die Mess- und Auswerteeinheit (18) 
einen Impedanz-Analysator (64) umfasst. . 

13. Vorrichtung nach einem der Anspruche 1 bis 12, dadurch 
gekennzeichnet, dass die Mess- und Auswerteeinheit (18) 
mit einer Mess- und Steuersof tware ausgestattet ist, 
die gewonnene und/oder abgeleitete Messdaten an eine 
relationale Datenbank ubergibt, die vorzugsweise mit 
einer Auswertesof tware verknupft ist. 

14. Vorrichtung nach Anspruch 13, dadurch gekennzeichnet, 
dass die Auswertesof tware eine Fitf unkt ionalitat zur 
Berechnung theoretischer Impedanzspektren fur die ein- 
zelnen Materialproben umfasst, wobei die Berechnung 
vorzugsweise unter Zugrundelegung eines Schaltkreis- 
aquivalents (90) erfolgt, das mindestens ein elektroni- 
sches Bauteil umfasst. 
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15. Vorrichtung nach Anspruch 13 Oder 14, dadurch gekenn- 
zeichnet, dass die Auswertesof tware eine Datamining- 
Funktionalitat umf asst - 

16. Vorrichtung nach Anspruch 15, dadurch gekennzeichnet , 
dass die Datamining-Funktionalitat unter Anwendung vor- 
zugsweise mehrdimensionaler Zielf unktionen arbeitet. 

17. Vorrichtung nach einem der Anspruche 14 bis 16, dadurch 
gekennzeichnet, dass die Datamining-Funktionalitat eine 
Visualisierungsf unktionalitat umf asst . 

18. Vorrichtung nach einem der Anspruche 1 bis 17, gekenn- 
zeichnet durch eine Heizeinrichtung (22), in welche die 
Probenplatte . (12) vorzugsweise eintauchbar ist. 

19. Verfahren zur Analyse einer Probenplatte (12), auf der 
mindestens zwei Materialproben angeordnet sind, umfas- 
send folgende Schritte: 

- Messen eines Impedanzspektrums fur jede der Material- 
proben; 

- Ablegen der gemessenen Impedanzspektren in einer Da- 
tei Oder einer Datenbank; 

- Bestimmen eines Aufbaus eines Schaltkreisaquivalent s 
in Abhangigkeit von dem jeweils gemessenen Impedanz- 
spektrum fur jede der Materialproben, wobei das jewei- 
lige Schaltkreisaquivalent mindestens ein elektroni- 
sches Bauteil, insbesondere mindestens einen Widerstand 
und/oder mindestens ein RC-Glied umfasst; 

- Bestimmen von fur eine Fehlerminimierungsrechnung er- 
forderlichen Startwerten fur die Bauelemente des jewei- 
ligen Schalt kreis Equivalent s ; 
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- Berechnen eines theoret ischen Impedanzspektrums fur 
mindestens eine der Mat erialproben mittels der Fehler- 
minimierungsrechnung unter Zugrundelegung des fur diese 
Materialprobe gemessenen Impedanzspektrums sowie der 
Startwerte fiir die Bauelemente des betreffenden 
Schalkreisaquivalents ; 

- Bestimmen von Fitwerten fur die Bauelemente des 
betreffenden Schalt kreisaquivalents ; 

- Bestimmen einer ValidierungsgroBe fur das errechnete, 
theoretische Impedanzspektrum; 

- Bestimmen einer Auswertegroiie durch Vergleich mindes- 
tens eines der Fitwerte fur die Bauelemente mit einem 
Ref erenzwert . 

20. Verfahren nach Anspruch 19, dadurch gekennzeichnet , 
dass eine Anzahl an in Reihe geschalteter RC-Glieder 
unter Berucksicht igung eines vorzugsweise vorwahlbaren 
Schwellenwerts bestimmt wird, wobei vorzugsweise maxi- 
mal vier RC-Glieder ausgewahlt werden. 

21. Verfahren nach Anspruch 19 Oder 20, dadurch gekenn- 
zeichnet, dass die Startwerte fur die Bauelemente eines 
ersten RC-Glieds des Schalt kreisaquivalents in Abhan- 
gigkeit von der maximal gemessenen, imaginaren Impedanz 
Z"_iyiAX ermittelt werden, wobei ein Startwiderstand 
R1_START und eine Start kapazitat C1_START berechnet 
werden . 

22. Verfahren nach einem der Anspruche 19 bis 21, dadurch 
gekennzeichnet , dass die Fehlerminimierungsrechnung 
durch Variation der Dimensionierung der einzelnen Bau- 
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elemente des Schalt kreisaquivalents um vorzugsweise 1 % 
durchgef uhrt wird. 

23. Verfahren nach einem der Anspriiche 19 bis 22, dadurch 
gekennzeichnet , dass bei der Fehlerminimierungsrechnung 
ein Fehler des theoretischen Impedanzspektrums durch 
Analyse der Differenz zu dem gemessenen Impedanzspekt- 
rum ermittelt wird. 

24. Verfahren nach einem der Anspriiche 19 bis 23, dadurch 
gekennzeichnet , dass fiir jede der Materialproben Impe- 
danzspektren unter ver schiedenen Pruf gasatmospharen und 
vorzugsweise bei verschiedenen Temperaturen gemessen 
werden . 

25. Verfahren nach einem der Anspruche 19 bis 24, dadurch 
gekennzeichnet, dass die AuswertegroiSe fur jede Materi- 
alprobe in eine Datenbank geschrieben wird und anhand 
der in der Datenbank abgelegten Auswertegrol^en ein Da- 
tamining durchgefuhrt wird. 

26. Verfahren nach Anspruch 25, dadurch gekennzeichnet, 
dass das Datamining mittels einer Zielfunktion durchge- 
fuhrt wird. 

27. Verfahren nach Anspruch 26, dadurch gekennzeichnet, 
dass das Datamining mittels einer visuellen Datamining- 
Funktionalitat durchgefuhrt wird. 

28. Verfahren nach einem der Anspruche 19 bis 27, dadurch 
gekennzeichnet, dass die gemessenen Impedanzspektren 
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mittels einer Kontrollf unkt ionalitat visuell liberpruft 
und/oder ausgewertet werden . 

29. Datenverarbeitungsanlage mit einem Programm zur Durch- 
fuhrung eines Verfahrens nach einem der Anspruche 19 
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